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1 Introduction 

Solar system bodies that lack a significant atmosphere and internal magnetic 
fields, such as the Moon and asteroids, can to a first approximation be considered 
passive absorbers of the solar wind. The solar wind ions and electrons directly 
impact the surface of these bodies, and a wake is created behind the object. In 
a frame that moves with the solar wind plasma, the refill of the lunar wake can 
be viewed as the more general problem of plasma expansion into a vacuum. 

One-dimensional plasma expansion into a vacuum has been studied for a 
long time, see [1] and references therein. The specific case of one-dimensional 
plasma expansion into the lunar wake has also been studied [2 [3] . A disadvan- 
tage of one-dimensional models is that the magnetic field component along the 
dimension is forced to be constant, from the requirement that V • B = 0. On 
the other hand, three dimensional models [I] does not allow the same spatial 
resolution and size of the simulation domain as lower dimensional models, due 
to computational constraints. Thus, a two-dimensional model of the lunar wake 
refill could be useful. 

A one-dimensional approximation of the lunar wake is an initial state with 
uniform plasma, corresponding to the solar wind, surrounding a region of vac- 
uum (zero density) of width equal to the lunar radius. With time the plasma 
will expand into the vacuum region. This approximates plasma conditions along 
a line that at t = lies in the terminator plane and goes through the center 
of the Moon. The line then convects with the solar wind as time goes by. The 
magnetic field at t = corresponds to the interplanetary magnetic field (IMF) 
in the solar wind. It is possible to use the same convecting frame approach also 
in two dimensions. If we consider the terminator plane at t = and let this 
plane convect with the solar wind plasma bulk velocity, u sw , it will at time t 
be at x = —tu sw . The initial configuration at t = is the yz-plane with proton 
number density n = uq for r < R, and n = n sw for r > R. Here the cc-axis 
is anti-parallel to the solar wind, and r = \J y 1 + z 2 . We study such plasma 
expansion into a vacuum in a two-dimensional geometry using a hybrid model 
with R = 1730 km. The hybrid model has ions as particles and electrons as 
a mass- less fluid, and is described in 0], and references therein. In the hybrid 
model the electric field computation involves a division by charge density. Thus, 
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there is a problem in vacuum regions. Here we handle this by letting no have 
a small, non-zero value. In this case uq — 0.1n sw . However, this low density 
population of ions is not shown in the plots. 

We can note that more complicated low dimensional models are possible. 
We could follow the convecting plane (or line in the one-dimensional case) from 
upstream, in the undisturbed solar wind, and after each time step remove ions 
that collide with the Lunar surface from the ion velocity space distribution. 
This removal would continue until the plane has convected past the Moon. The 
solution would however no more be independent on the solar wind velocity. 



2 Results 

In what follows, the solar wind number density, n sw = 5 cm~ 3 , the ion tem- 
perature is 5 • 10 4 K, the electron temperature is 10 • 10 4 K, and the IMF is in 
the xy-p\an.e with a magnitude of 5 nT. For these plasma parameters, the ion 
inertial length is 100 km, the Alfven velocity is 49 km/s, the thermal proton 
gyro radius is 60 km (the gyro time is 13 s), and the ion plasma beta is 0.35. 
The cell size is 78 km and the number of meta particles per cell is 170. 

How the wake number density depends on the direction of the initial mag- 
netic field is shown in Fig. [I] The time t = 200 s of the plot corresponds to a 
distance of tu sw downstream of the moon, e.g., for u sw = 500 km/s the distance 
is 10,000 km w 5.8 Lunar radii. 

We see a density rarefaction that expands from the Moon at the fast mag- 
netosonic speed, as has been seen before in three-dimensional simulations [4] . It 
does not propagate parallel to the IMF as can be seen for the case of an IMF at 
a 90 deg angle to the solar wind flow direction. At this time, for these condition, 
the wake has partially refilled except for the case of a deg IMF (parallel to the 
solar wind flow) for which we still have a (deformed) vacuum region. We also 
see that the wake is flattened in the plane of the IMF due to the fact that the 
ions flow along the IMF. 

For the 90 deg IMF case we see density oscillations in the refilled region, 
where surrounding plasma intrudes. In the next section we derive criteria for an 
instability in the hybrid model that can be a possible cause for these oscillations. 
Note that these oscillations are a feature of the hybrid model, it is unknown if 
such oscillations also would occur in a 2D particle in cell model of the Lunar 
wake that includes also electrons as particle. 

3 Ion beam instability in a hybrid model 

For longitudal waves propagating along the magnetic field lines, the hybrid 
model reads 

^ + U M + ^M = (i) 

dt dx rrii dv 
where the electric field is obtained from the electron momentum equation as 

E= MkJL(»y /3 , (2) 
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Figure 1: Proton number density in the 1/2-plane at t — 200 s for an IMF direction 
of 90, 45, 20, and degrees away from the image normal (solar wind direction). The 
plane of the IMF is horizontal. The black circle show the outline of the initial low 
density region of radius 1730 km (representing the Moon). 
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and quasi-neutrality gives 

n = J ,fd Z v (3) 
Linearization = /o + /i , E = E\, n = n + rii gives 

ot ax mi 
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and 

ni = / /id 3 f (6) 



Laplace transformation <9/<9i -4- — zcj and 9/9x -4- ifc gives 



+ + = o, (7) 
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and 

= • • /7,-/ M - ifc / / ;( /Y = //,•/.•, / ' ' : / v. 
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(8) 

Eliminating the common factor E\ and rearranging the terms, we obtain the 
dispersion relation 

1 SW| i &t 

3 rriin J v — oj/k ov 

Here the frequency and wavenumber appear only as a ratio to/k. If we replace 
this ratio by v p , the "phase velocity", we obtain the equation 

3 TTLiUQ J V — Vp ov 

Note that in this equation, v p depends only on the plasma parameters but not 
on lu or k. Hence v p is a complex number, v p = vr + ivi, which can either have 
positive or negative imaginary parts vi. The dispersion relation can be written 
lu = (vr + ivj)k. If vj is positive, then we have an instability which formally 
goes to infinity as k — »■ oo. In a simulation, however, the wavenumber is limited 
by the grid size. 
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